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Bose- Einstein correlations of identically charged pion pairs were measured by the PHENIX experi- 
ment at mid-rapidity in Au-|-Au collisions at ^snn = 200 GeV. The Bertsch-Pratt radius parameters 
were determined as a function of the transverse momentum of the pair and as a function of the cen- 
trality of the collision. Using the full Coulomb correction, the ratio flout /Aside is smaller than unity 
for (fcx) from 0.25 to 1.2 GeV/c and for all measured centralities. However, using recently developed 
partial Coulomb correction methods, we find that -Rout/-Rsido is 0.8-1.1 for the measured (/st) range, 
and approximately constant at unity with the number of participants. 

PACS numbers: 25.75.Dw 

Following its origin in the study of proton-antiproton extensively used to measure source distributions in rela- 
annihilations P| , Bose-Einstein correlations have been tivistic heavy ion collisions 0, |^ • These measurements 
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were originally motivated by theoretical predictions of 
a large source size and/or a long duration of particle 
emission 0, 0, 0| ~ which would result from a soften- 
ing of the equation-of-state in a first-order phase transi- 
tion to a quark-gluon plasma (QGP). The technique of 
Bose-Einstein correlations is based upon quantum sta- 
tistical interference, but final state interactions such as 
Coulomb repulsion modify the relative momentum distri- 
butions for pairs of identical particles emanating from the 
collision region. Both effects are included in multidimen- 
sional Gaussian fits to the normalized relative momentum 
distributions yielding the fit parameters which are the 
RMS-widths in each dimension, i?iong, -Rside, -Rout 
also referred to as HBT radii. A finite duration of emis- 
sion leads to an effective elongation of one of the HBT 
radii, i?out, which is parallel to the mean transverse mo- 
mentum of pair, fcx = (pit + P2t)/2; however, the other 
radii along the axes perpendicular to fcx represent the 
geometrical source sizes. Therefore the i?out/-Rsidc ra- 
tio could become larger than unity even for a cylindrical 
particle source, if a long duration of particle emissions 
occurs. For dynamic (i.e. expanding) sources, the HBT 
radii depend on the fcx, and correspond to lengths of ho- 
mogeneity, regions of the source which emit particles of 
similar momentum . Measuring the fcx dependence of 
the HBT radii provides essential constraints on dynam- 
ical models that include the space-time evolution of the 
source 

Hydrodynamic models for the space-time evolution of 
a rehadronizing QGP predicted that the measurement of 
the i?out/-Rsido ratio at moderate values of fcx provides 
a sensitive test of the long duration of particle emission, 
a signal of a slowly burning first order phase transition 
from QGP to hadrons at RHIC 0. Their predictions 
that the i?out/-Rsidc ratio should reach ~1.5 at of 
^^0.5 GeV/c were not borne out by initial measurements 
of Bose-Einstein correlations of pions from Au-|-Au col- 
lisions at ^snn = 130 GeV This disagreement 
between theory and data has been called the "RHIC HBT 
puzzle" H. 

We present here data on Bose-Einstein correlations of 
charged pion pairs measured by the PHENIX experiment 
at RHIC for Au-|-Au coUisions at ^snn = 200 GeV. In 
this analysis, we adopt a recent fitting technique that 
provides for a self-consistent treatment of the Coulomb 
final state interaction for a source that is made up of a 
smaller core and a more extended halo of long-lived res- 
onances. We introduce a new parameterization in which 
the strength of the Coulomb interaction is constrained 
by the measured unlike-signed pion correlation. 

The PHENIX detector provides particle identification 
(PID) capabilities for hadrons, leptons and photons over 
a wide momentum range. The setup of the PHENIX de- 
tector has been described in detail elsewhere . In this 
analysis, we use the west arm of the central spectrom- 
eter, which covers the pseudorapidity region \rj\ < 0.35 



and A(j) = tt/2 in azimuthal angle over a region of 0.2 
GeV/c < fcx < 2.0 GeV/c. The drift chamber (DC), 
at a radial distance between 2.0 m and 2.4 m, provides 
trajectory information in the azimuthal direction. A pad 
chamber (PCI) at 2.5 m provides z-coordinate informa- 
tion. Combining the DC and PCI information, a track 
model provides a 3-dimensional trajectory and momen- 
tum vector for charged particles. The momentum res- 
olution is Sp/p ~ 0.7% ® 1.0% X p (GeV/c), where the 
first term is due to the multiple scattering before the 
DC and the second term is the angular resolution of the 
DC. For this analysis, the electromagnetic calorimeter 
(EMCal) provides the time of arrival of particles at its 
front face located 5.1 m from the beam axis. The timing 
resolution is approximately 400 psec for hadrons. This 
analysis is based on a sample of 34 million minimum 
bias events taken with a magnetic field of 0.78 T-m and 
triggered by the coincidence of the Beam-Beam Coun- 
ters (BBC) and Zero-Degree Calorimeters (ZDC) - cor- 
responding to 92 ± 2% of the total inelastic cross section 
of 6.8 b. Event centrality is determined from the corre- 
lation between the BBC multiplicity and the analog re- 
sponse of the ZDC. About 23 million events are selected 
with a requirement that the collision vertex measured by 
the BBC has \z\ < 30 cm. Each track is required to have 
an associated hit on the EMCal within 2<t of the track's 
projection to the EMCal, where cr refers to the resolu- 
tion of the projection. Charged particles are identified 
by the time-of-fiight technique using timing information 
between the BBC and the EMCal, combined with mo- 
mentum and flight length calculated by the track model. 
Charged particles in the PID zone within 1.5 a of the 
ideal squared-mass peak of pions but 1.5 a away from 
the kaon bands are identified as pions. After the track 
quality and PID cuts, ~45 million positive pions and ~51 
million negative pions are selected in a momentum range 
from 0.2 to 2.0 GeV/c. 

The pion correlation function is experimentally defined 
as C2(q) = A(q)/_B(q), where A{q) is the measured two- 
pion (actual pair) distribution of pair momentum differ- 
ence q, and B(q) is the background pair (mixed pair) 
distribution generated using mixed events from the same 
data sample. Event mixing is done selecting events that 
have similar multiplicities and event vertices. As in our 
previous analyses at ^snn = 130 GeV actual pairs 
within 1 cm in the beam direction (AZoc) and 0.06 ra- 
dians in azimuthal angle {A(j)Dc) in the DC are elimi- 
nated from the pair sample to remove ghost tracks, then 
within 5 cm in AZdc and 0.03 radians in A(f)DC are also 
eliminated to avoid a strongly inefficient region for pairs. 
Since pairs which are close to one another in the EMCal 
are either reconstructed as single hit clusters or are af- 
fected by cluster-sharing, pairs whose hits are within 8 
cm at the EMCal are eliminated. The background mixed 
pairs are subject to the same cuts as actual pairs. Af- 
ter the pair cuts, ~110 million positive pion pairs and 
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~140 million negative pion pairs remain. The number 
of pairs is 40 times larger than the data s amp le used for 
the PHENIX data at = 130 GeV 0. The sys- 

tematic errors on the HBT radius parameters from the 
pair cuts are evaluated by varying the pair cuts to be 
~2% for i?side and i?iong, and ~7% for i?out- The correc- 
tion for the multi-track reconstruction efficiencies in the 
DC and EMCal are determined by a GEANT-based ^ 
Monte Carlo (MC) simulation of the detector. The pair 
efficiencies of the DC and EMCal are estimated by the 
ratio of the distribution of separations of actual pairs to 
that of normalized mixed pairs. Also the multiplicity de- 
pendent pair efficiencies are estimated by an embedding 
technique: simulated pion pairs are embedded in real 
events and the reconstruction efficiency for these simu- 
lated pion pairs versus multiplicity is calculated. The 
pair reconstruction inefficiencies in the DC and EMCal 
are corrected by using pair efficiency factors estimated by 
the embedded MC simulation. Finally, we removed pairs 
within 0.005 radians in A(j)oc where a pair inefficiency 
still remains even after the correction by MC. Correction 
for the residual HBT effect jT^ is estimated as a system- 
atic error. The acceptance for this analysis is large, as a 
consequence the systematic error is ~1% on each radius 
parameter. 

In order to compare directly to previous analyses, we 
apply the standard "fuW strength Coulomb correction 
calculated from the pair Coulomb wave function ^ for 
a 3-D Gaussian parameterization of the source using 
the radii obtained from the previous fit. We fit the 
Bose-Einstein correlation with the full Coulomb correc- 
tion to the the 1-D qinv parameterization, C2((/inv) = 
1 + Ainv exp(— _Rj?jjyq?jy), and the 3-D Bertsch-Pratt pa- 
rameterization is given by 

C2 = 1 + Aexp(-i?^.^pq3^jj^ - -Rout^out - -Rlong^long)- (1) 

The relative momentum q is decomposed into 
(Qsidc, (Zout, 9iong), whcrc the longitudinal component 
(qiong) is parallel to the beam axis, the out component 
(^out) is parallel to the mean transverse momentum of 
the pair, and the side component (gside) is perpen- 
dicular to both gioiig and (/out0j 111- This analysis is 
performed in the Longitudinal Center-of-Mass System 
(LCMS), where the mean longitudinal momentum of 
the pair vanishes. In this frame, the duration of particle 
emission couples exclusively to gout- In a general case, 
cross-terms may appear in Eq. ^ but they vanish in our 
measurement of central collisions at mid-rapidity due to 
symmetry reasons |l8| . 

In the realistic source, however, many of relatively 
long-lived particles (e.g. 77, rj') which decay into pions 
and have a Bose-Einstein interference too narrow to be 
resolved by experiment also have a Coulomb interaction 
that is negligible. For this reason, a new functional form 
was proposed |0| which is based upon a Core-Halo pic- 
ture of the source l20| and assumes that the fraction of 



pairs. A, which have Bose-Einstein interference are the 
only pairs that contribute to the Coulomb interaction. 

C2 = Ccore + Chalo = [A(l + G)F] + [1 - A] (2) 

where G corresponds to the Gaussian term in Eq. ^ 
and F is the Coulomb correction term. This correction 
method was recently adopted by CERES 21] and STAR 
[2^ . However, our formula differs slightly, because the 
momentum resolution correction to A (which is small ac- 
cording to MC simulation) is included in our systematic 
error, rather than being incorporated into F. 

In order to test the underlying hypothesis of Eq. |21 

we fit the strength of the Coulomb interaction (A^ ) 

to the unlike-signed correlation function in the range 
0.2 < /ct < 2.0 GeV/c, and obtained a value for 
A+_ = 0.50 ± 0.04 with x^/DoF = 3.0. This value is 
clearly inconsistent with the full strength Coulomb cor- 
rection but greater than the value of A, especially for 
fcx < 0.5 GeV/c. We attribute this difference mainly 
to the uj resonance, which is sufficiently long-lived to be 
unresolved with Bose-Einstein correlations, but may con- 
tribute significantly to the Coulomb interaction. To ac- 
count for this contribution, we have modified Eq. [21 to 

fix the total coulomb strength at A_| while allowing the 

Bose-Einstein fraction. A, to vary from zero to A^ , 

C2 = [A(l + G)F] + [(A+_ - A)F] + [1 - A+_] (3) 
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FIG. 1: Panels (a) and (b) show one-dimensional correla- 
tion functions for tt^tt^ and 7r~7r~. The bottom figures show 
the three-dimensional correlation function for 7r~7r~ with the 
full Coulomb (opened circle) and without Coulomb (filled tri- 
angle) corrections for 0.2 GeV/c < k-r < 2.0 CeV/c for 0- 
30% centrality. The projection of the 3-D correlation func- 
tions are averaged over the lowest 40 MeV in the orthogonal 
directions. The error bars are statistical only. The lines over- 
laid on the open circles (filled triangles) correspond to fits to 
Eq. (Eq. |5J over the entire distribution. Panel (c) shows 
the one-dimensional correlation function of unlike-signed pi- 
ons for 0.2 < fcT < 2.0 GeV/c. The two overlaid histograms 
show calculations for the full (dashed) and the 50% partial 
(solid) Coulomb corrections. 
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In applying this formula, we still calculate the addi- 
tional Coulomb fraction using the Bertsch-Pratt source 
of approximately 5 fm, rather than estimating the larger 
source distribution for the a; decay products. Therefore 
we use this formula to provide an upper bound on the 
effect of the additional Coulomb interaction. The differ- 
ence between Eq. O and Eq. |2 is used in our estimate of 
the systematic errors. 

Figure ^ shows the one-dimensional correlation func- 
tion for 7r+ pairs, pairs, and unlike-signed pion pairs 
along with projections of the three-dimensional correla- 
tion functions onto gsidc, ^out and giong for 7r~ pairs. 

Figure 121 shows the fcx-dependence of A, Bertsch-Pratt 
radii, and the ratio i?out/^?sido for the 30% most cen- 
tral events, corresponding to (A'part) = 281. For the 
full Coulomb correction, A is approximately constant in 
all the fcx bins while the radius parameters fall rapidly 
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FIG. 2: The At dependence of the Bertsch-Pratt radius pa- 
rameters and A for charged pions for 0-30% centrahty. Filled 
triangles show the results from fits to a core-halo structure by 
Eq. 121 with statistical error bars and systematic error bands. 
Open circles and squares show the results from the full (Eq. 
and 50% partial (Eq. j^J Coulomb corrections with statisti- 
cal error bars, respectively. Results at 130 GeV by PHENIX 
are given by filled circles. 



with increasing fcx- The full Coulomb corrected radius 
parameters at ^snn = 200 GeV are slightly different 
from the full Coulomb corrected radius parameters at 130 
GeV 01 at the same (A;T)because of the improved pair 
efficiency correction, but similar within errors. For the 
partial Coulomb correction, results from the two differ- 
ent correction methods show similar trends. Compared 
to the full Coulomb correction i?side and i?iong systemat- 
ically decrease while i?out increases, and A decreases in 
the low fcx region. In the case of the full Coulomb correc- 
tion, the ratio i?out/-Rside, in Fig. [3 is around 0.6-0.8 up 
to fcT'~^l-2 GeV/c. On the other hand, i?out/-Rsido from 
the partial Coulomb correction is systematically larger 
than that from the full Coulomb correction, and slightly 
decreases from ~1.1 to ~0.8 as fcx increases. Reanalysis 
of the 130 GeV/c data with Eq. (21 gave results that are 
fully consistent with the 200 GcV/c results. 

Figure shows the collision centrality dependence 
of the radius parameters. The number of participants 
[Npart) is evaluated from the charged particle multiplicity 
using a Glauber model calculation [23j. To evaluate the 
-^part dependencies of the Bertsch-Pratt radius parame- 
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FIG. 3: Bertsch-Pratt radius parameters versus the cube root 
of the number of participants (A^p^ft) for charged pions for the 
fits to a core-halo structure by Eq. [51 (filled triangles) with 
statistical error bars and systematic error bands, in the 0.2 
GeV/c < fcT < 2.0 GeV/c range, with (/ct) = 0.45 GeV/c. 

1/3 

The dashed lines show fits to po -|- Pi * -^part- Fitted po and 
pi values with fits to the core-halo stracture are given in each 
panel. The opened circles show results with the full Coulomb 
correction by Eq. Q 
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1 /3 

ters, we fit with a function of po + Pi * ^part- Foi" the 
/wH Coulomb correction, the fit pi parameters indicate 
that i?side and -Riong show similar iVpart dependencies and 
i?out has a slightly smaller iVpart dependence. For the fits 
to a core-halo structure, all radius parameters show sim- 
ilar A^part dependencies. All radii are consistent with a 

1 /3 

linear increase with A^p^rt • The i?out / -Rsidc ratios are ap- 
proximately constant for all centralities. The ratios from 
the partial Coulomb corrections are systematically higher 
than those using the full Coulomb corrections. These fit 
parameters are given in Fig. 21 

In conclusion, we have presented the Bertsch-Pratt 
HBT radii in the LCMS for identified charged pions mea- 
sured by PHENIX in Au-|-Au collisions at ^/snn = 200 
GeV. The fcx dependence of the HBT radii was mea- 
sured for (A'^part) — 281, and the centrality dependence 
was measured for At = 0.45 GeV/c. These measure- 
ments are consistent with results from Au-f Au collisions 
at 130 GeV when a similar analysis {full Coulomb cor- 
rection) is performed. We also performed two differ- 
ent partial Coulomb analyses, one based upon a self- 
consistent treatment of the Coulomb correction, and the 
other based upon direct comparison to the unlike-signed 
correlation, which is shown to be inconsistent with the 
application of a full Coulomb correction. Both Coulomb 
corrections (Eq. [21 and Eq. (SJ yield similar values of 
-Rout/^sido which slightly decreases from ~1.1 to ~0.8 in 
the range of fcx = 0.2 - 1.2 GeV/c for (iVpart) = 281, 
and approximately constant at unity with the number of 
participants for (/ct) = 0.45 GeV/c. These detailed mea- 
surements of the transverse momentum dependence of 
the HBT radii, in particular that of i?out/^side, provide 
extremely strong constraints for model builders. 
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